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The elementary excitations of strongly correlated one-dimensional electronic systems - spinons and holons - are discussed in 



o 



T3 

O 



in 
o 



i 

T3 



an exact solution of the Haldane-Shastry and Kuramoto-Yokoyama model. We derive and exactly solve the equation of 
motion both for the two-spinon wavefunction and for the one-spinon, one-holon wavefunction. By solving the equations of 
motion we find the spinon-spinon and spinon-holon interaction to be a short-range attraction. The physical consequences of 
such an attraction on the spin- form factor and on the hole-spectral function are also worked out. 
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I. INTRODUCTION 



Phenomenological Landau's Fermi liquid theory [1] describes many of the condensed matter systems. The basic 
assumption of Landau's theory is that the spectrum of the interacting Fermi liquid may be adiabatically deformed 
^ ' to the spectrum of a noninteracting Fermi gas. Equivalently, one may say that the interaction strength can be 
continuously increased from (Fermi gas) to its Fermi Liquid value without hitting a singular point. The excitations 
I ' of a Fermi liquid can be continuously deformed to electrons and holes and are referred to as quasiparticles and 
, quasiholes, respectively. 

In a "non-Fermi liquid" such a scenario does not apply. Typical examples of non-Fermi liquids are the one- 
dimensional strongly correlated electronic systems. The elementary excitations of these systems are collective modes 
carrying charge 1, but no spin ("holons"), and modes carrying spin- 1/2, but no charge ("spinons"). Originally found 
as elementary excitations of spin-1/2 one-dimensional antiferromagnetic spin chains [2,3], spinons appear as localized 
£f~) ' spin-1/2 spin defects, embedded in an otherwise featureless spin-singlet spin-liquid state. Spinons and holons carry 
"fractions" of the quantum numbers of spin waves and/or holes. 

In these lectures we analyze the dynamics of spinons and holons in one-dimensional electron-systems close to half 
filling. In particular, we study their interaction in the framework of simple, exactly solvable model, using a formalism 
different from the usual Bethe-ansatz approach [4] . Such an approach treats collective modes of strongly correlated 
systems (spinons, holons, etc.) as quantum mechanical particles. In particular, it is possible to associate a two-body 
wavefunction to a two-spinon or to a spinon-holon pair. The corresponding equation of motion can be straightforwardly 
solved. The form of the interaction potential in the "Schrodinger" operator and the behavior of the solution as a 
O ■ function of the separation among the particles indicates a short-range attraction. 

^ , In the thermodynamic limit, the spinon-spinon and spinon-holon attraction drives the low-energy physics of one- 
dimensional antiferromagnets and strongly-correlated insulators. It is responsible for the sharp features at threshold 
in the spin structure factor, and in the hole spectral function. These features have been experimentally observed [5,6] 
' and they provide an important evidence for a short-range attractive force. 
The notes are organized as follows: 

• In Section II we introduce the Haldane-Shastry model of a one-dimensional antiferromagnet. We study in some 
details its basic features, its ground state and its elementary excitations; 

• In Section III we discuss the two-spinon wavefunction, together with the corresponding equation of motion, its 
solution, and physical interpretation. We prove the existence of a short-range attraction among spinons; 

• In Section IV we consider one-holon and one-spinon states, within the framework of the Kuramoto-Yokoyama 
model which is the supersymmetric extension of the Haldane-Shastry model; 
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• In Section V we study the one-spinon one-holon wavefunction. We solve the corresponding equation of motion, 
and study in detail the basic features of the solution. Finally, we prove the existence of a short-range spinon-holon 
attraction; 

• In Section VI we spell out the physical consequences of the attraction between spinons and holons. We trderive 
its effects on the dynamical spin structure factor in one dimensional antiferromagnets and on the angle-resolved 
hole spectral function in one-dimensional insulators; 

• In Section VII we provide our conclusions. 



II. THE HALDANE-SHASTRY MODEL OF A ONE-DIMENSIONAL ANTIFERROMAGNET AND ITS 

ELEMENTARY EXCITATIONS 

Spinons appear as elementary excitations of one-dimensional half-odd spin antiferromagnets. We explore the dy- 
namics of spinons in the framework of the Haldane-Shastry (HS) model of a spin-1/2 antiferromagnetic chain [7,8]. 
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FIG. 1. a) Haldane-Shastry interaction; b) Spin-spin static correlations in the ground state of the Haldane-Shastry model. 

The HS model is a system of spin 1/2 particles lying at the sites of a lattice with periodic boundary conditions. 
The lattice may be thought of as "wrapped" around a circle. The interaction among the spins is antiferromagnetic, 
with interaction strength J, and inversely proportional to the square of the chord between the corresponding sites. 
Let N be the number of lattice sites, so that a generic point on the lattice will be labeled by a = 1, . . . , N. Let us 



also define the analytic coordinate on the lattice, z a = exp 



2ma/N 



The HS-modcl Hamiltonian reads: 



J/2tA 2 " 

9 I M I 2^ 



Za - Z f3 \ 7 



(1) 



where S Q is the spin at site a ((S Q ) 2 = 3/4). 

Due to the magic of analytic variables on the circle, z* = l/z a , the interaction in Eq.(l) is, in fact, an analytic 
function of the z a : s. This allows us to work out the action of Hhs on many-body states, depending on the z a 's, in 
terms of differential operators acting on the analytic extension of the corresponding wavefunctions. 

Hhs clearly has the global S'C/(2)-symmetry generated by the total spin S = X^a=i ^ a ^ so P osses an additional 
symmetry, corresponding to total spin-current conservation and generated by the operator A, given by: 



JV 



l_ \ ~> / Zg + Zfj 

2 2-? \z a -Z8 



(S Q x S^) 



(2) 



that is, one has: 



[Whs.S] = [Hhs, A] = 



The action of Has on the ground state and on one- and two-spinon collective states of the system will be the subject 
of the remaining part of this Section. 
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A. Ground state of Tins 



Let |^gs) be the ground state of Tins- Due to the low-dimensionality, in ID quantum antifcrromagnets with short- 
range correlations, we cannot expect |\&gs) to order. Rather than being realized as a one-dimensional antifcrromagnctic 
state with Neel order, the ground state of the system is given by a one-dimensional spin-liquid state, with spin density 
uniform over length scales of the order of the lattice step. The ground state also has short-range correlations and no 
net magnetization. For the HS model, in the even- N -case, I^gs) is expressed in terms of its components on states of 
the Fock space with N/2 spin-|, the remaining ones being J,. In particular, let \z±, . . . , z N / 2 ) be defined as: 

N_ 

\z u ...,z„) = f[S+\l,...,l) (3) 
i=i 

where | j, . . . , J.) is the fully polarized states with all the spins at the lattice sites being J., while z\, . . . , z N / 2 are 
iV/2-lattice coordinates. We define |\I/gs) by: 

J A 

2 2 

(z 1 ,...,zn\^ gs ) = ^ gs (z 1 ,...,zn)= Y[ (z>- Zj ) 2 Y[z t (4) 

i<j = l t=l 

Tins can be traded for a differential operator, when acting on functions of the z Q 's as discussed in detail, for instance, 
in [9]. Here we outline the main idea and briefly sketch the basic steps of the derivation. 

First of all, we have to split the scalar product S a • S/3 into (S+Sp + S~Sp)/2 + S^Sp and to separately consider the 
action of the various operators on the ground state. Therefore, let us observe that [S^ S^^gs](zi, . . . , Zm) is different 
from zero only if one of the arguments z\, . . . , zjv/2 equals z a . By Taylor expanding the corresponding function, we 
obtain: 

(where {zj} is a shorthand for {zi, . . . , Zm})- 

The coefficients A\ are calculated in [9,10]. In particular, it is possible to show that they are all zero for N > I > 2. 
Therefore, we can rewrite Eq.(5) as: 

™f (N-l)(N-5) _ _N-3 2 d_ 1 3 0M f 9os(z hlll2 ZM)\ 
£1 12 3 2 Z >dz 3 + 2 Z idz]j{ z 3 } 



where, in deriving Eq.(6), we have used the identity: 



( N M 11 

T^kF^r s(zi "-' ZM) • (6) 

v n-Ah ' J ' y 



V 2 Z 2 7 2 



(7) 



(z a - zp){z a - Zj) (zp - z a )(z0 - Zj) (z~, - z a ){z 1 - Zp) 
Finally, it is straightforward to prove that: 

Eqs.(6,7) lead us to the final result: 

H HS ^gs{zu---,zn ) = E gs ^gs(zi,...,zn) (8) 

with: 
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that shows that ^gs{zi, ■ ■ ■ ,zm) is, in fact, an eigenstate of Tins- By employing an approach similar to the one 
leading to Eq.(8), it is also possible to show that ^gs is a total spin singlet: [S^gsK^i, • ■ • i z m) — 0. The proof that 
\1/gs is the only ground state of 7Yhs is not immediate. It involves the "factorization property" of 7Yhs , first discovered 
by Shastry [11], and discussed at length in Ref. [10]. We will not prove it here, but will rather focus on some relevant 
physical properties of the state described by \I>gs- 

In Fig. lb), we report the spin-spin correlations in ^gs, (^Gs|Sa ■ S^I^gs), as a function of the separation, a — (3. 
We may clearly recognize the algebraic decay of spin correlations which is expected for a spin- 1/2 chain where the 
spectrum of the elementary excitations on top of the ground state shows no gap [12]. As the total spin momentum 
of the state is zero and the spin density is uniform we can refer to "Jgs as to a disordered spin-singlet, homogeneous 
nondegenerate spin liquid. Elementary excitations on top of such a spin-liquid spin-singlet are provided by collective 
modes, carrying total spin 1/2, the spinons. The next subsection deals with these particles. 



B. One- and Two-spinon states of the Haldane-Shastry model 

Spinons are stable excitations of the spin- 1/2 HS model. Therefore, they maintain their integrity upon interaction. 
This means that, in a scattering process, the numbers of incoming and outgoing spinons are equal. Therefore, it 
is possible to diagonalize Hhs within a subspace at fixed spinon number. As the total spin is also a constant of 
motion of the model, subspaces with fixed total spin polarization of the spinons do not mix one with each other upon 
interaction. This justifies the procedure we make use of in this Section, diagonalizing Hhs in the one-spinon and in 
the two-spinon fully polarized subspaces, respectively. 

In order to construct one-spinon states, let us consider the HS-chain with an odd number of sites. In this case, 
the minimum possible spin for each state is be 1/2. In particular, the ground state of the system will now be a 
superposition of states with total spin 1/2. The relevant spin- 1/2 states may be constructed from uniform spin-liquid 
states, by holding the spin at a site a to be J., for instance. As a function of the t spin-coordinates, z\,...,Zm 
(M = (N — l)/2) , the corresponding wavefunction reads: 

M M M 

k p (zi, ...,z M )= - *<*) n (* - *i? n * w 

3=1 i<j=l t=l 

In Fig. 2a), we sketch the spin density and the charge density profile for the state in Eq.(9). As it appears from the 
plot, 'I'qP describes a uniform spin- liquid state, just like \&gs does, except for a localized "spin bump" at z a , carrying 
total spin -1/2. Such a spin bump is a "localized J,-spinon at z a " . 

Localized one-spinon states are not translationally invariant. Therefore, they cannot be eigenstates of the Haldane- 
Shastry Hamiltonian. Eigenstates of Has ma Y be constructed, instead, by propagating the spinon, that is, by using 
functions in the form: 

/v 

1 AT _ 1 

^(z u ...,z M ) = -J2«r^(zi,...,z M ) ;0<m<^- (10) 

a=l 

Indeed, by following the same steps leading to Eq.(8), it is possible to show that the eigenvalue equation for \I/® P 
becomes: 

|*»p = A* S Q P (11) 
^From Eqs. (10,11), we derive the eigenvalue E m , corresponding to 



N-3 



M - z D 



dz n 
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Egs + 2" 



51 -(« 



sp)2 



where, in Eq.(12), we have introduced the crystal momentum = ttN/2 — 2-n{m + 1/A)/N 
Let us discuss some important features of the one-spinon dispersion relation, E(q^) = 



_ J 
2 



sp\2 



(12) 



First 



of all, let us notice that the one-spinon Brillouinc zone (BZ) is halved, that is, — n/2 < < n/2. This corresponds 
to the absence of negative energy one-spinon states and, correspondingly, to the halving of the number of available 
states. Secondly, as we may see also from the plot of the dispersion relation reported in Fig. 3a), the dispersion relation 
becomes gapless at the endpoints of the BZ. The existence of gapless excitations is a typical feature of half-odd spin 
antifcrromagnets [12], and is the reason for the algebraic decay of the spin-spin correlations we see in Fig. lb). 
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a) 



b) 



FIG. 2. a) Profile of the spin- and charge-distribution of a localized spinon at z a \ b) Profile of the spin- and 
charge-distribution of a localized holon at ho- 




a) b) 

FIG. 3. a) One-spinon dispersion relation; b) Negative-energy part of the one-holon dispersion relation. 

Since spinons maintain their integrity, we can build multi-spinon states by first creating M-localized spinons and 
then by propagating each one of them. A residual spinon-interaction will require a further diagonalization of the 
Hamiltonian matrix restricted to the multi-spinon plane-wave subspace, as we will see in the particular case of two- 
spinon states. 

To construct the state for two localized spinon at z a and zp, let us consider a HS-chain with even N . Let the 
M(= N/2 — 1) spins at Z\,...,Zm be | and let the spin at z a and zp be held 1. The corresponding collective 
wavefunction will be given by: 

M MM 

...,z M )= H(zj - z a )(z 3 - z fi ) n - zjf n Zt (13) 

j=l i<j=l t=l 

A set of independent two-propagating-spinon wavefunctions is given by the states: 
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JV N 

*% n (z u ...,z M ) = -S ^E(W*> ^) ;m>n = O...M (14) 



AT2 

a=l 8=1 



The Schrodinger equation for \I>^ now reads: 



7? 



7V-3 



2 

In the plane-wave basis, Eq.(15) becomes: 



d d 
z <*x h z /3"5 — 

OZ a OZ0 



+ [2M 2 -M(iV-2)]|* s Q p /3 = A^ p /3 (15) 



J/2tt\ f AT 2 / 19 24 



2 V -/V / 1 48 V AT 



+m(M - to) + n(M n) -^JL _ £ (m _ n + 2*)*£ +/>n _ / (16) 

1 > i=i 

where Im = n if n + m < M, Im = M — to, if n + m> M. Eq.(16) has solutions of the form [13,10]: 

Im 



=0 



where the coefficients a™" are recursively defined as: 



I O 

m - n + a ™« = l (18) 

2^ + m-n+i)^ 



The corresponding eigenvalues are given by: 

Emn = E GS + 



J /2tt x : 
2\N 



fur \ \ m — n 

m(M - to) + n(M - n) - J — - — - 



(19) 



Eq.(19) shows that the energy of the two spinons is equal to the ground state energy plus the energies of each single 
spinon plus an interaction contribution that is subleading in the thermodynamic limit. Indeed, as A^ — > oo, we obtain: 

E mn - E GS -+ E(q%,q°v) = E(q%) + E(q%>) (20) 

The fact that, in the thermodynamic limit, the energy is additive does not mean that the effects of spinon interac- 
tion disappear as A~ — > oo. Instead, as we are going to see in the next Section, spinon interaction safely survives 
the thermodynamic limit, and appears to be the driving mechanism of the low-energy physics of one-dimensional 
antiferromagnets. 

III. SPINON DYNAMICS IN THE HALDANE-SHASTRY MODEL 

In this Section we study spinon interaction and its properties. To do so, we have to introduce a formalism that 
allows to treat spinons as actual quantum mechanical particles. 

In Section II, we referred to the states as the two-spinon energy eigenstates, and to the states as the 
states for two localized spinons at z a and z@. Therefore, it is always possible to fully decompose any of the ^^'s 
within the set of the 's: 



ran 

M 



*«0 = E E *a OTn(?)*£n (21) 
m— n— 



G 



where V!^ n (z) is a polynomial in z of degree m — n. 

By definition, the coefficients of the combination in Eq.(21) provide the coordinate representation for the wavefunc- 
tion for two spinons in a state of energy E mn . The corresponding equation of motion for V^ n can be derived from 
the following identity chain (which is a direct consequence of Eq.(15)): 

E mn {^ n \^%) = «JWhsTO> 



(M 



_d_ _d_ 

OZ a OZ a 



+(M - z —)z 9 1 ( Za + Zf} 

Zf3 dzf3 Zl3 dzp 2 \z a - zp 



--<><" 



' dz, 



mn I a/3 



(22) 



In the differential operator at the r.h.s. of Eq.(22) we recognize the energy operators of the two spinons, plus an 
interaction term that is velocity dependent and diverges as 1/x at short separation x between the two particles. 

From the equation ($^J*^,) = {®*& n \® a £ n ) z™ zJP% n {^) and from Eq.(22), we may derive the differential equation 
for the relative coordinate wavefunction: 



1 

m 

2 



n — (— — m + n)z 



d_ 

dz 



r^ n (z) = o 



(23) 



Eq.(23) is a degenerate hypergeometric equation, with parameters a = — m + n, b = 1/2, c = a + b. Its solution is 
provided by the hypergeometric polynomial [14]: 



T~'mn( z ) — T^m-nyZ) 



r[m - n + 1]_ ^ T[k + i]r[m 



r[i]r[m-n + i] ^ r[fc + i]r 



fc + i] 



21 z k 



(24) 



l^mnC 6 * 1 )! 2 gi yes the density of probability for the configuration with two spinons at a distance x, as a function of 
x. In Fig. 4 we plot [P^ n {e tx )\ 2 vs. x. While at large separation the probability density oscillates and averages to 1, 
as it is appropriate for noninteracting particles, when the two spinons come close to each other, a huge probability 
enhancement appears utill the function reaches its maximum for x — 0. The physical meaning of such a behavior may 
be summarized in two basic points: 

• Spinon interaction is short ranged: This is clearly seen when the probability amplitude approaches the one 
for two noninteracting particles in the case when the spinons are far from each other; 

• Spinon interaction is attractive: This is seen in the huge probability enhancement for configurations when 
the spinons are on top of each other. 

Therefore, we have identified spinon interaction as a short range attraction. We now turn to the question of the fate 
of such an interaction in the thermodynamic limit. In order to do so, let us look at Fig. 4b), where we show an insert 
of the plot of \P^ n {e lx )\ 2 around x = 0, for the same relative momentum, but for increasing N. We see that, the 
larger N is, the stronger the enhancement. Therefore, not only spinon interaction safely survives the thermodynamic 
limit but its effects become more relevant as the the sample size grows. 



60 







iPVh 



71 



500 



-0.1 






a) b) 
FIG. 4. a) Square modulus of the two-spinon wavefunction for N = 100 and m ~ n ~ N/2 — 1 as a function of spinon 
separation; b) Inset, at zero spinon separation, of the square modulus of the two-spinon wavefunction for increasing N. From 
bottom to top: N = 200, 400, 600, respectively. 
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IV. THE KURAMOTO-YOKOYAMA MODEL OF A ONE-DIMENSIONAL INSULATOR AND ITS 

ELEMENTARY EXCITATIONS 



In this Section we discuss the Kuramoto-Yokoyama (KY) model for a strongly correlated one-dimensional electron 
system. The KY-model is the supersymmetric extension of the HS-model. Strong on-site Coulomb repulsion forbids 
double occupancy, but one may have fewer electrons than sites. This allows for charged excitations (holes). 

The Kuramoto-Yokoyama Hamiltonian is defined as [15]: 

Hky = k (I) r^f Pg { s " ' 8/3 " 5 ? &* c *> + \^+^)-^- i} p o 

where c aa is the single electron annihilation operator, n a = J2a c aa c aa is the charge operator, and the Gutzwillcr 
projector: 

N 

pg = n i 1 - °U c ii c ci c ^) 

a=l 

accounts for the no-double occupancy constraint. 

In addition to the usual "bosonic" symmetry, generated by the total spin operator S and by the total charge 

operator N — ^2 a n a , Hky also possesses four Fermionic symmetry operators, given by Q a — X^ = i c a(7 ano - by 

ei = Ea=i4. [9]. 

Exactly at half- filling, Hky reduces back to Has- Therefore, the ground state of Hky at half- filling is still given 
by I^gs)- Also, charge-0, spin-1/2 spinon excitations will be described by the same wavefunctions \P^ P and we 
have derived for one spinon in the HS-model, with the same dispersion relation. 



A. One-holon and one-holon one-spinon states of the KY-model 

In order to construct charge-1, spin-0 holons, we may consider a state where one electron is removed from the top 
of a spinon, so that the resulting charge vacancy will carry no spin. The dispersion relation of such an excitation 
spans, in general, the whole Brillouin zone. Therefore, both positive- and negative-energy holons may arise, as charged 
excitations of the KY-model. Because of the supersymmetry of Hns > a positive-energy holon state may be obtained 
by acting with Qj on one-spinon fully polarized states. As those states will not be relevant for our analysis, we will 
neglect them in the following, and will rather focus onto wavefunctions with negative-energy holons. 

In the case of odd N, the wavefunction for a propagating negative energy holon is given by: 

N N N /V + 1 

. . . , z M \h) = h n H( Zj - h) J] (* - \{zf^<n< —I- (25) 

j=l i<j=l t=l 

where M = (N — l)/2, and z\, ... , zm denote the position of the |-spin electrons, h is the position of the empty site. 

Clearly, carries total charge 1. By using the technique introduced for the HS-model, it is straightforward to 
prove that S Zl t>^ = S~^° = 0. Therefore, S z ^° is a spin-0, charge-1 state, that is, a one-holon spin singlet. 

The mathematical derivation of the action of Hky on is long and boring, although straightforward. Here, we 
will not go through the various mathematical steps, but will rather quote the final result: 

[Hky*!HW,...,;»,I'>> = {-./QQ ( w "^)+^(f) »(n-^±i)}*5°(*„...,*„|/.) 

M- J (fi)( W + ^^) +£( ' r)) }*"° <26) 

In Eq.(26) we have introduced the one-holon crystal momentum, = irN/2 + 2n(n — 1/A)/N (mod 27r), and the 
dispersion relation: 

(27) 



E(ql°) 
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The negative-energy part of the dispersion relation in Eq.(27) is plotted in Fig.3(b). Again, one may recognize that 
the dispersion relation becomes gapless at the corners of the Brillouin zone, similarly to what happens for spinons. 
The remaining half-part of the Brillouine zone is spanned by positive-energy states, which we are not considering 
here. 

The state for a localized holon at ho, provided by the Fourier transform of ^>^° reads: 



^ho \p /,-™\T/ h ° 



(28) 



n=0 



Let us now introduce the one-spinon one-holon states of the KY-model. The starting point is provided by the states 
*n?s h °7 where the holon is propagating, while the spinon is localized at s. Let N be even, M = N/2 — 1, and let: 



M 



M 



M 



h n n ife - s )fe - h )\ n & - z o? n * ; i < « < m + 2 



j = l i>j=l 

The KY-Hamiltonian acts on *^ ho as follows [16,17]: 



-E GS ]^°( Zl ,...,z M ) = 

2 V N ' 



(29) 



n-M) + (M-l)4- S 292 ' Ik + Sf 9 



ds ds 2 2 



+ 



hs 



{h-sf 



,T,Sp,ho _ ( £\" 1 ^Sp^O 
\h) ^ S = h 



(30) 



In order to construct the corresponding energy eigenstates of Hkm, let us consider the one-spinon, one-holon plane 
waves: 



*^ ho (zi, ...,z M )=J2 ^^(zi, ...,zm) 



(31) 



In the plane wave basis, the eigenvalue equation takes the form: 



2 V N ' 



sp,ho 
ran 



N 

m(M — to) + n(n — 1 + — ) 



sp,ho 



if m — n + 1 < 0, and: 



1 m 

--(n - m + l)*^ h ° - (n - m - 1) 2 



Hky 

2 V N I 



m(M — to) + n(n — 1 



^ mn 



ho 



(32) 



M—m 



.hm-n + lW^°-{m-n+l) £ 



if to — n + 1 > 0. 

Eqs. (32,33) take the following eigenvectors and eigenvalues: 



(33) 
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• Case to — n + 1 < 0: 



sp,ho 



sp,ho 
m-j,n+j 



E + = E GS + -(— 



AT 1 

m(M — m) + n(n — 1 — — ) — -(n — m — 1) 



(34) 



Case to — n + 1 > 0: 



M—m 



r[i + |] 



\p / ,y x LV "T- 2 J ,T,sp,ho 

^ l J r[i]rb' + i] 



^mn "(jo 2 V iV 



AT 1 

to(M — to) + n(n — 1 — — ) + — (n — to — 1) 



(35) 



Notice that, in either case, the energy can be written in terms of the spinon and holon momenta as: 

E mn = E GS + E(q%) + E(q*°) - <£°| 



(36) 



From Eq.(36) we see that, as for the case of two spinons, in the thermodynamic limit the energy of a spinon-holon 
pair on top of the ground state energy reduces to the sum of the energies of the single isolated spinon and of the single 
isolated holon. As for spinons, this does not imply the interaction between spinons and holons has no effects as the 
size of the system goes large. 

V. ONE-SPINON ONE-HOLON DYNAMICS IN THE KURAMOTO-YOKOYAMA MODEL 

In order to write down the Schrodingcr equation for a spinon-holon pair, Eq.(23), in analogy with two spinons, we 
need to define the state of a localized spinon at s and a localized anyon at ho- Fourier transforming 1 i' l n p ^ back to 
coordinate space performs this task: 



M+2 



sp.ho 

,ho 



( Zl ,...,z M )=Yl h o n K P . 



ho 



(37) 



71 = 1 



sp.ho 

s,h 



can be fully decomposed within the basis of one-spinon, one-holon energy eigenstates. In particular, we obtain: 



■ s,h 



M+2 n-2 M+2 M 

E E sm v i ^ ho (^«„ ho +E E * m ^o n ^ (£» ] 



n—l m—0 



n—1 m—n — 1 



(38) 



Eq.(38) defines the one-spinon one-holon relative wavefunction V^ n ho - To write down the equation of motion for the 
one-spinon one-holon wavefunction, we need the string of identities: 



2 V N > 



-ran IV j ( 2n \2 



k ° 1 __ } /$Sp,h0|^,Sp hox 

7 I \^mn I * s,h I 



s- h \h 



(39) 
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where v — M \i m — n + 1 < 0, u = otherwise. 

From Eq.(39), it is straightforward to derive the differential equation for P^ n ho : 



1 



2±- 

dz (1 - z) 



y M-m+l 



(I-Z) 



r s ^ ho (z) = o 



if m — n + 1 < 0, and: 



-T5 S P 
' ran 



ho 



(l) r 



V^ n ho (z) = 



(40) 



if m — n + 1 > 0. 

Notice that, unlike the spinon case, Eqs.(40) are first-order differential equations. The reason for this is that 
spinon and holon energy bands have opposite curvature and, because of supersymmetry, the second derivatives of the 
dispersion relation have equal absolute values, but opposite signs. Nevertheless, one-spinon one-holon dynamics is 
basically the same as two-spinon dynamics, as we are going to discuss next. 

Let us solve Eqs.(40). It is easy to see that they allow for polynomial- like solutions, given by: 



"P S P 
' ran 



ho 



(*) = 



M—m—1 

E • 

fc=0 



r[fc- 



r[i]r[* + i] 



z k (m - n + 1 < 0) 



and: 



T[k+i] .1 



-(-) fe 



(m- n + 1 > 0). 



(41) 



The squared modulus of '(e* x ) gives the probability of finding a spinon and a holon at a distance x from each 
other. Such a probability is plotted in Fig. 5a). Notice that, as for a spinon pair, at large separations between the two 
particles, the probability oscillates and averages to 1 flatly. This is a signal of the absence of long-range interaction 
effects. As the two particles get close one to each other, on the other hand, the probability peaks up, and shows a 
consistent enhancement at zero separation. 

Therefore an interaction between spinons and holons exists. Such an interaction is short-ranged and, since it favors 
configurations with the two particles on top of each other, it is attractive. In conclusion, we assert that, as it happens 
among spinons, there is a short-range attraction as a relevant interaction between spinons and holons. The physical 
consequences of such an attraction is the subject of the next Section. 



iPIV)i 



500 




CV)i 2 



8000 




-0.1 
b) 

FIG. 5. a) Square modulus of the one-spinon one-holon wavefunction for iV = 100 and m — n ~ N/2 — 1 as a function of 
spinon-holon separation; b) Inset at zero spinon-holon separation of the square modulus of the one-spinon one-holon wavefunc- 
tion for increasing N. From bottom to top: iV = 100, 200, 400, respectively. 



VI. PHYSICAL CONSEQUENCES OF THE INTERACTION BETWEEN SPINONS AND HOLONS 

In this Section we address the question of how spinon-spinon and spinon-holon interaction can be seen in a real 
experiment. 

The experiments we think about are neutron scattering measurements performed on quasi one-dimensional antifer- 
romagnets [5], and ARPES spectra measurements, performed on quasi one-dimensional insulators [6]. In the former 
case, one measures the dynamical spin form factor, that is, the imaginary part of the spin susceptibility, defined as: 
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A spin (u,q) = l%m{(y GS \S(-u; 7 -q) ■ S(w, g )|M> G S>} 

In the latter case, on the other hand, the photoemission of an electron from the sample is induced and, by detecting 
the emitted photoelectron, one eventually reconstructs the spectral function of the recoiling photohole, defined as: 

A holc (uj,q) = -^m{G holc (^,q)} 

where Ghoie(w,g) is the hole Green function. 

By using the results we derived in the previous Section, we can exactly calculate the contribution to *4 spin (w, q) from 
two spinon states, as well as the contribution to *4 hole (w, q) coming from one-spinon one-holon states, A sp ' ho (oj, q). In 
particular, we can express these contributions in terms of the two-spinon and of the one-spinon one-holon wavefunc- 
tions, calculated when the two particles are on top of each other. In this way, we directly see what the measurable 
consequences of the probability enhancement (and therefore the interaction between particles) are. 

The wavefunction for propagating spin-one excitation with momentum q = 2nk/N is given by: 

N m 

S^ GS = $>;) fe * s a p « = JV £ £(-ir + "< i+ "<5(n + m k)V sp n (l)^ p n (42) 

a m— n— 

Therefore, the form factor will be: 

\ ^ ran I ran! ^ 



{iV lit 
EE 
m— n n=n 



From Eq.(43) we see that A sp (u>, q) is fully determined only by the squared modulus of the two-spinon wavefunctions at 
zero separation between the two spinons, that is, by the probability enhancement. Since the probability enhancement 
is a direct consequence of spinon attraction, we get to the ultimate result that the spin form factor in one-dimensional 
antiferromagnets is fully determined only by spinon attraction. 

An expression similar to Eq.(43) holds for A sp ' ho (uj,q). In particular, we obtain [16,17]: 

rM+2 1-2 



a (^-^4^2^ + 10+ {E + n _ Egs) 



where, again, we see that the one-spinon one-holon contribution to the hole spectral function is entirely determined 
only by the probability enhancement, that is, by spinon-holon attraction. 

The relevant quantities appearing in Eqs. (43,44) are the squared modulus of the wavefunctions when the two 
particles are on top of each other, and the norms of the many-body two-spinon and one-spinon one-holon states. The 
values of those quantities are reported in Appendix A, where we show that they are basically given by products and 
ratios of Euler's T- functions. Therefore, taking the thermodynamic limit of Eqs. (43,44) is just matter of properly 
using Stirling's approximation: 

T[z] w V^(z - lj'-ie-**- 1 ) (z large) 
As discussed in detail in Refs. [10,17], the final result is: 

A spin {uj q) = J @ ^ 2 ^ ~ ~ u; +i(g)] Q [ a; ~ (45) 

4 ^[u - Zn^Xqj\^J\u) - u+i{q)\ 

(where tO-i{q) — J/2q(ir — q), uJ+i(q) = J/2(2ir — q)(q — ir), and ^2(3) = J/2q(2n — q)), for the spin form factor, and: 



12 



2 Jq\j u- J\q-\\ 2 



1 /J[g+f] 2 -c 



(46) 



Tr- 



ior the one-spinon, one-holon contribution to the hole spectral density. 

The quantities reported in Eqs. (45,46) share some important features. Indeed, we see that no resonances appear 
in either case. A resonance would be an evidence for a stable spin-one propagating spin-wave in the former case, 
and for a stable propagating hole state in the latter case. We rather see broad features at the tails of the spectra, 
basically signaling the nonintegrity of the spin-wave versus decay into a spinon pair, and of the hole versus decay 
into a spinon-holon pair. Moreover, the spectra take sharp spikes at the creation threshold, that is, when the two- 
spinon/spinon-holon pair is created at the minimum of energy. Such a sharp threshold appears quite a natural 
consequence of our results, that is, it is the effect of the attraction, taken in the thermodynamic limit. 

Sharp threshold features have been seen in the experiments quoted above [5,6]. They should be interpreted as a 
direct experimental evidence of attraction among spinons and spinons and holons. Such an attraction appears to be 
ubiquitous for fractionalized excitations in one-dimensional strongly correlated systems, despite the fact that we have 
derived it within the framework of a somehow oversimplified mathematical model. 



In these notes we have reviewed some results concerning the interaction among collective excitations of strongly cor- 
related one-dimensional electronic systems: holons and spinons. A detailed mathematical derivation, performed within 
the framework of the Haldane-Shastry and of the Kuramoto-Yokoyama model, the simplest models of one-dimensional 
antiferromagnet and of one-dimensional strongly correlated insulator, provided us with important informations about 
such an interaction and its nature. We have seen that two spinons, as well as a spinon and a holon, interact by means 
of a short-range attraction. Although such an attraction creates a probability enhancement that goes large in the 
thermodynamic limit, it is not able to create a two-spinon or a one spinon-one holon bound state (a propagating 
spin wave or a propagating hole, respectively). The effect of the attraction is rather the creation of a sharp emission 
threshold in the spin-wave and in the hole spectral function, on top of a broad two-spinon, or one-spinon one-holon 
continuum. Such a sharp threshold is seen in neutron scattering experiment on quasi one-dimensional antiferromag- 
nets [5], and in ARPES spectra measurements on quasi-one-dimensional insulator. In our view, it provides relevant 
evidence for an attractive interaction between spinons, and between spinons and holons. 



In this Appendix we provide the formulas for the norms of two-spinon and one-spinon, one-holon states and 
for the two-spinon and one-spinon, one-holon probability enhancements. We will skip the long and boring, though 
straightforward, mathematical derivation of the results of this Appendix, and refer the interested reader to the original 
papers [10,17]. 

The starting point for the derivation of {&mn\^mn) ^ s a formula derived by K. Wilson [18]: 



VII. CONCLUSION 



APPENDIX A: NORMS AND PROBABILITY ENHANCEMENTS 




(Al) 



The next step is to consider the recursion relations 



($£nl$£n) _ (n-\){M-n+^)(m-n+lf 



(*m-i,J*m-i,n> <M - n + l)(m - n + §)(m - n + \) 



and: 




/<J> s p |<J> s p \ 

\ mn I mnl 



(M — m + l)(m — n + \){m — n — 




(A2) 



— II m,n — II 



(M - m + \){m + \){m - n) 2 



Eqs.(A2) are solved by: 



13 



^ SP i^p v = r r[m - n + |]£[m - n+ j] £[m + 1]T[M - m+ |] 
\ mn\ van) r 2 [m-n+l] r[m + |]r[M - m + 1] 



I>+ |]r[M-n+l] 



r[n + l]r[M - n + §] 



(A3) 



where C M may be calculated from Wilson's integral in Eq.(Al), and is given by C M = N M (2M)\(M + ±)/(tt2 m ). 
A similar technique may be applied, to derive (&mn°\^mn°)- The result is: 

($S p,Ho |$S p >ho) = N M + Ami (M + i r[M-m + i]r[m + i] 



for m — n + 1 < 0, and: 

<*r°l*r°> = iV- +1 »(M + l) r j M - m+ y r [m + 3 1 ! (A4) 

^ V 2 ; T[M-m+ l]r[m+ |] 



for m — n + 1 > 0. 

The probability enhancement V^ n (l) may be derived from standard properties of hypergeometric functions [14]. 
The result is: 

rfiir[m-n + il 

1 [m — n + 2J 

Finally, V^^°(l) may be recursively derived, as discussed in the Appendix of Rcf. [17], and the result is: 

, , rfM-m + i] 
^ ho (i) = 2-L- 



r[i]r[M-m] 

for m — n + 1 < 0, and: 



v\m + -1 



for m — n + 1 > 0. 
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